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The ability to control electric currents in solids using diodes and transistors is undoubtedly at the origin of 
the main developments in modern electronics which have revolutionized the daily life in the second half of 
20th century. Surprisingly, until the year 2000 no thermal counterpart for such a control had been proposed. 
Since then, based on pioneering works on the control of phononic heat currents new devices were proposed 
which allow for the control of heat fluxes carried by photons rather than phonons or electrons. The goal of 
the present paper is to summarize the main advances achieved recently in the field of thermal energy control 
with photons. 


I. INTRODUCTION 

The diode and the transistor introduced by F. BraurP 
and Bardeen et al^ (based on the works of Julius E. 
Lilienfeld from 1925), in 1874 and 1948 are the main 
building blocks of almost all modern systems of infor¬ 
mation treatment. These elementary devices allow for 
rectifying, switching, modulating and even amplifying 
the electric current in solids. Similar devices which 
would provide the same degree of control of heat cur¬ 
rents instead of the control of electric currents are not 
as widespread in our daily life and the concepts for such 
devices for heat flow management at the nanoscale were 
introduced as recently as about ten years ago. In 2006 
Baowen Li et alW have proposed a thermal counterpart 
of a field-effect transistor. In this device the tempera¬ 
ture bias plays the role of the voltage bias and the heat 
currents carried by phonons play the role of the electric 
currents. Later, several prototypes of phononic thermal 
logic gateJ^ as well as thermal memorie^^ have been de¬ 
veloped in order to process information by phononic heat 
currents rather than by electric currents. 

However, this technology suffers from some weaknesses 
of fundamental nature which intrinsically limit its per¬ 
formance. One of the main limitations comes proba¬ 
bly from the speed of acoustic phonons (heat carriers) 
which is four or five orders of magnitude smaller than 
the speed of light and therefore also some orders of mag¬ 
nitude smaller than the speed of electrons. Another in¬ 
trinsic limitation of phononic devices is related to the 
inevitable presence of local Kapitza resistances. These 
resistances which originate from the mismatch of vibra¬ 
tional modes at the interface of different elements can 
reduce the phononic heat flow dramatically. Finally, the 
strong nonlinear phonon-phonon interaction mechanism 
makes the phononic devices difficult to deal with in pres¬ 
ence of a strong thermal gradient. On the contrary the 
physics of energy transport mediated by photons instead 
of phonons remains unchanged close and far from ther- 
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mal equilibrium. These limitations and difficulties ex¬ 
plain, in part, why so many efforts have been deployed, 
during the last decades, to develop full optical or at least 
opto-electronic architectures for processing and manag¬ 
ing information. Particularly important developments 
have been carried out during the last decade with plas- 
monics system^^H ^ith the goal to increase the speed of 
information processing substantially while reducing the 
dimension of devices to the nanoscale at the same time. 

In this paper we review the recent developments made 
to achieve a control of heat flow with contactless devices 
using thermal photons. We discuss the operating princi¬ 
ples of the three basic building blocks for heat flow man¬ 
agement which are the photonic thermal diod^, pho¬ 
tonic thermal transistoi^ and photonic thermal mem- 
or based on phase-change materials. 

II. RADIATIVE THERMAL DIODE 

Asymmetry of heat transport with respect to the sign 
of the temperature gradient between two poi n^ts Fig , [^a) 
is the basic definition of thermal rectificatiorP^^ which 
is at the heart of a variety of applications as, for ex¬ 
ample, in thermal regulation of buildings as scetched in 
Fig.^b). Usually, this dissymmetry in the thermal be¬ 
haviour of a system finds its origin in the dependence 
of the electical/phononic/optical properties of materials 
with respect to the temperature. The effectiveness of the 
thermal rectification is generally measured by means of 
the normalized rectification coefficient which can be de¬ 
fined as 

r, = 

max(4>^, 4>^) 

where and denote the heat flux in the forward 
and backward operating mode, respectively. Different 
solid-state thermal rectifiers have been suggested d uring 
the last decade based on nonlinear atomic vibrationJi^EIl, 
a nonlinear dispersion relation of the electron gas in 
metalJI^, direction dependent Kapitza resistance^^^ or 
based on the dependence of the superconducting density 
of states and phase dependence of heat currents owing 
through Josephson junction^^. 
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FIG. 1. (a) Sketch of thermal rectification principle. Forward 
scenario (top-left) : the left side of a medium is in contact 
with a hot reservoir at temperature Tl while its right side 
is in contact with a cold reservoir at temperature Tr < Tl 
so that a heat flux is exchanged between both reservoirs 
through the medium. Backward scenario (top-right) : the 
temperature gradient is reversed resulting in a heat flux 
in the opposite direction. Due to a strongly nonlinear thermal 
behavior of the medium the magnitude of flux and for 
can be different resulting in a rectification of the heat flux in 
one preferred direction, (b) Application of radiative heat flux 
rectification to a passive thermal regulation of a building. The 
non-linear behavior of insulating materials allows to reduce 
both heating/cooling during summer and winter. 


During the last five years, photon-mediated thermal 
rectifiers have been proposed both in planai^I^Jt^ and 
non-planar geometr}^^ By several authors to tune ra¬ 
diative heat exchanges both in near field (i.e. for sep¬ 
aration distances smaller than the thermal wavelength) 
and far field (i.e. for separation distances larger than 
the thermal wavelength) using materials with tempera¬ 
ture dependent optical and scattering properties. How¬ 
ever, in plan geometry, only relatively weak thermal rec¬ 
tifications coefficients were ac hieve d with the proposed 
mechanisms (77 < 44% in RefsP^^, 77 < 52% in Ref.l^. 
Unfortunately, these rectification coefficients are not suf¬ 
ficiently high to operate such devices as radiative diodes. 
Very recently another rectification mechanism driven by 
the phase-change of material properties was proposed 
in order to increase significantly the rectification coef¬ 
ficient both in near-fielcP^ and far-held regim^. In 
such insulator-metal transition (IMT) materials, a small 
change of the temperature around its critical temperature 
Tc causes a sudden qualitative and quantitative change 
of the material properties from an insulating to a metal¬ 
lic behaviour by a Mott transitiorP^. Thanks to the IMT 
the optical propertieP^of such materials undergo a rapid 
change as well resulting in thermal rectihcation coeffi¬ 
cients as large as 60% in far-held regime (Fig. | 2 }b) and 
99.9% in near-held regime (Fig. [^c). 

To illustrate this, let us consider a system as depicted 


in Fig. where two plane samples, one made of VO 2 
and one made of amorphous glass (Si 02 ) at tempera¬ 
tures Tl and T^^, respectively. Both media are sepa¬ 
rated by a vacuum gap of thickness d. Here, VO 2 is a 
IMT material which has a critical temperature of about 
Tc = 340 K. For temperatures below Tc VO 2 is in an 
insulating phase and has a monocilinic structure and for 
temperatures above Tc it is in a metallic phase with a ru¬ 
tile structure. The optical propertied can be modelled 
by a uni-axial permittivity for T < Tc and by a isotropic 
permittivity for T > T^. 

We now examine this system in the two following ther¬ 
mal operating modes: (i) In the forward mode (F) the 
temperature Tl = Tc+AT of the VO 2 slab is greater than 
its critical temperature Tc so that VO 2 is in its metallic 
phase while the temperature of the glass medium is set 
to Tr = Tc — AT < Tl. Obviously the average temper¬ 
ature is then centered around Tc and the net radiative 
heat flow is from the left-hand side to the right, (ii) In 
the backward mode (B) we choose Tl = Tc — AT so that 
VO 2 is in its insulating phase and Tr = Tc + AT > Tl- 
Again the average temperature is centered around Tc but 
the net energy flow is this time from the right-hand side 
to the left. 

The exchanged radiative heat flux per unit surface be¬ 
tween t wo semi- infinite media can be written in the gen¬ 
eral fornP^J^^^ 


^F/B-J U 7 7^:5^FF/B(w,K;d) 

° i={s,p} 

POO 

= / dcj A0 (cj)(/:?f/b(<^5 

Jo 


( 2 ) 


where A0(ci;) = 0(ci;,Tl) — 0(co’,Tr) is the difference of 
mean energies of Planck oscillators at frequency uj 


0(u;,T) 


1 


exp 


huj 

kBTi/j 


)-i) 


(3) 


at the temperatures of two interacting media; /cr is Boltz¬ 
mann’s constant and 2Tvh is Planck’s constant. As for 
represents the energy transmission coeffi¬ 
cient of each mode (cj, n) {k, is the wave vector parallel 
to the interfaces) for the two polarization states (s and 
p polarization). For in general anisotropic m edia i t is 
defined for the propagating modes {k < uj/c) 


TJ'T/b = Tr[(l - IR][^IRr)Dlr(1 - HlHl)]Dlr^] (4) 


and for the evanescent modes {k > uj/c) hy 

7',f/b = Tr[(]Rj, - Rr)Dlr(1EIl - 

(5) 

introducing the component of the wave vector in the vac¬ 
uum gap normal to the interfaces ^ . Note 

that the difference between 7 },f and is caused by 
the choice of the temperatures for the two modi (i) and 
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(ii). The reflection matrices of each interface depend on 
these temperatures through the optical properties of the 
media. They are given by (/ = L, R) 




rr(c.,K) ry 


The matrix is defined as 


©lr = (1 - ^ 


(6) 

(7) 


The matrix elements of the reflection matrix are the 
reflection coefficients for the scattering of an incoming j- 
polarized plane wave into an outgoing j'-polarized wave. 
For isotropic or uniaxial media where the optical axis 
is orthogonal to the surface there is no depolarisation, 
i.e. we have = 0. The remaining reflection 

coefficients are then given by 


kzO ^l\s 

^zO T ^l]s 


( 8 ) 


^p,p ^ ^i,\\kz0 ki-,p .g. 

^z,|| ^zO T ki-,p 

where ki-s^p are solutions of the Fresnel equatiorP^ 

“ ^*.11 = 0 - 

( 10 ) 

Here and ei^± are the permittivities parallel and per¬ 
pendicular to the surface of the uniaxial material VO 2 
assuming that the optical axis of VO 2 is normal to the 
surface. For amorphous glass which is isotropic we have 
ey = = esi02 • 

With the above expressions we can determine the heat 
flux in the forward and backward case for different tem¬ 
perature differences AT. The resulting rectification co¬ 
efficients T] of the system in Fig. [^a) are plotted in Fig¬ 
ure Hb) in far-held regime. It decreases monotonously 
with the temperature difference AT between both sam¬ 
ples. For very small temperature differences we see that 
T] ^ 60% illustrating the high rectihcation efficiency of 
the phase transition of VO 2 . This efficiency is very large 
given the fact that we use the simplest possible geome¬ 
tric. This prediction has been recently experimentally 
verihecP^. This rectihcation process could probably be 
father improved by texturing the medium as suggested 
in Ref.l^. 

In order to understand this large diherence of the ra¬ 
diative heat how in forward and backward direction we 
have plotted in Fig.j^the spectral heat hux (p{uj, d 00 ) 
dehned in Eq. (§ using the transmission coefficient Q 
for the forward and backward scenario. We note that 
when VO 2 is in its metallic phase (forward scenario, 
Tl > Tc), the spectral heat hux is broadband and scales 
approximately like oc . The same is true when VO 2 
is in its insulating phase. The diherence is that in the 


V02-a Si02 (a) 



Tl=T,+AT Tr=T,-AT Tl=T,-AT Tr=T,+AT 



FIG. 2 . (a) Phase-change radiative thermal diode. In the for¬ 
ward scenario (top-left) a IMT phase-change material (here 
VO2) is in its metallic state (V02-a) at a temperature higher 
than Tc. In the backward scenario (top right), the phase- 
change material is in its insulating state (VO2-C). (b) Recti¬ 
hcation coefficient in far-held regime of a V02/Si02 system 
with respect to the temperature diherence AT. The phase 
transition of VO2 occurs at Tc = 300 K. (c) Heat transfer 
coefficent h evaluated at T = Tc for VO2 is in its metallic or 
insulating state, respectively, h is plotted with respect to the 
separation distance d. 


insulating phase the spectral heat hux is generally larger 
and shows much more structure. This can be understood 
by the rehectivity of VO 2 which is plotted in Fig. |^b). 
In the metallic phase the rehectivity of VO 2 is relatively 
high as can be expected from a metallic material whereas 
in the insulating phase the rehectivity is generally much 
smaller in the shown frequency range. Therefore it fol¬ 
lows directly from Kirchhoh’s law that the emissivity of 
VO 2 in the metallic phase is smaller than in the insulat¬ 
ing phase. This asymetry is the key for obtaining a highly 
efficient phase-change thermal diode. It is interesting to 
note, that in the spectral region around 10^^ rad/s the 
rehectivity of VO 2 in its insulating phase can also be rel¬ 
atively large and attain values close or even larger than 
that in its metallic phase. In this frequency band, the 
insulating VO 2 is metal-like which means here that the 
permittivity can be negative. In this regime VO 2 in its 
insulating phase supports surface polaritons. 

In near-held regime that is to say when the separation 
distance between the VO 2 and Si 02 samples is smaller 
than the thermal wavelength, the ’diodicity’ of system 
becomes even better as we see on Figure |^c) by compar¬ 
ing the heat transfer coefficients 


hp/B 


Sffip/B 

dT 


( 11 ) 


for the forward and backward scenarios as a function of 
the separation distance d when temperature is equal to 
Tc. We observe a very large thermal contrast of about two 
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FIG. 3 . (a) Plot of the spectral heat flux (Pf/b (in far-held 
regime) introduced in Eq. § as a function of frequency for 
forward direction where VO2 is in its metallic phase and the 
reverse situation where VO2 is in its insulating phase. In 

2 

addition, we have plotted p — | for the case that both 

materials are perfect black bodies, i.e. 7 j,R/B = 1 for all prop¬ 
agating modes, (b) Plot of the rehectivity R of VO2 in its in¬ 
sulating and metallic phase for normal incidence as function 
of frequency. 


orders of magnitude in the near held against about a fac¬ 
tor 5 in the far held. In order to understand the physics 
involved in this spectacular increasement, we have plot¬ 
ted in Fig. the the transmission coefficients for both 
phases of VO 2 in the {uj^k) plane at a given separation 
distance in near-held regime. For insulating VO 2 we see 
in Fig.|^a) that this coefficient is close to one far beyond 
the light line k. = ujjc. In this region the contribution of 
the surface phonon polaritons of VO 2 which couple to the 
surface polaritons of Si 02 can be observed (see^for a de¬ 
tailed discussion). As the coupling is rather efficient in 
this region, the heat transfer is large. On the other hand, 
when VO 2 is metallic we obtain the transmission coeffi¬ 
cient shown in Fig. ib), the surface phonon-polariton 
contribution disappears and as a consequence the heat 
transfer is small. A detailed discussion of the efficiency 
of rectification in such diode as a function of the thick¬ 
ness of the slabs can be found in RefP. Also, based 
on a similar idea, a contactless thermal switch has been 
recently suggestecP^. 


FIG. 4 . (a) Transmission coefficient 7^,r/b (p polarization) 
between two insulating VO2 surfaces and (b) between two 
metallic VO2 surfaces. Both graphs were obtained for a sepa¬ 
ration distance d = 500 nm. The spectral character of thermal 
transfer changes radically across the phase transition 


III. PHOTONIC THERMAL TRANSISTOR 

The thermal diode allows the heat flow to propagate 
preferentially in one direction, but such a device does not 
allow for switching, modulating or even amplifying the 
flux exchanged between a hot and a cold body. How¬ 
ever, in recent years, physicists have begun to design 
and test thermal transistors with some succes^. In this 
paragraph we discuss the possibility to design a radia¬ 
tive thermal analog of electronic transistors as depicted 
in Fig. [^a). The classical field-effect transistor (FET) 
which is composed by three basic elements, the drain, 
the source, and the gate, is basically used to control the 
flux of electrons (the current) exchanged in the channel 
between the drain and the source by changing the voltage 
bias applied on the gate. The physical diameter of this 
channel is fixed, but its effective electrical diameter can 
be varied by the application of a voltage on the gate. A 
small change in this voltage can cause a large variation 
in the current from the source to the drain. 

Recently, a thermal counterpart of the FET has been 
proposed^ to control the near-field radiative heat ex¬ 
changes between two bodies. This near-field thermal 
transistor (NFTT) is depicted in Fig. [^b) and basically 
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consists of a source and the drain, labeled by the indices 
S and D, which are maintained at temperatures Ts and 
Td using thermostats where Ts > Td so that we have 
a net heat flux from the source towards the drain. A 
thin layer of IMT material labeled by G of width S is 
placed between the source and the drain at a distance d 
from both media and functions as the gate. This IMT 
material is able, as we have seen before, to qualitatively 
and quantitatively change its optical properties through a 
small change of its temperature around a critical temper¬ 
ature Tc. As in the previous paragraph we describe here 
the operating modes of radiative transistor using VO 2 as 
such IMT material. The choice of IMT depends on the 
operating temperature of the transistor. If the transis¬ 
tor should operate around T = 500 K then VO 2 could be 
replaced by LaCoOs, for instance. For the source and 
the drain we use again Si 02 . Hence we have in principle 
two oppositely biased diodes as discussed in the previous 
section which are connected in series so that the NFTT 
(for which we had a field-effect transistor in mind) can 
also be regarded as a bipolar transistor. 

Without external excitation, the system will reach its 
steady state for which the net flux received by the in¬ 
termediate medium, the gate, is zero by heating or cool¬ 
ing the gate until it reaches its steady state or equilibrium 
temperature T^. In this case the gate temperature Tq 
is set by the temperature of the surrounding media, i.e. 
the drain and the source. When a certain amount of heat 
is added to or removed from the gate for example by ap¬ 
plying a voltage difference through a couple of electrodes 
as illustrated in Fig. ib) or by extracting heat using 
Peltier elements, its temperature can be either increased 
or reduced around its equilibrium temperature Tq . This 
external action on the gate allows to tailor the heat flux 
between the source and the gate and the heat flux 
between the gate and the drain. 

In order to see, if this system can be operated as a 
thermal transistor we need first to determine the radia¬ 
tive heat flux in this system, which is a little bit more 
complex then the system studied for the radiative ther¬ 
mal diode. In a three-body system the radiative flux 
received by the drain takes the fornP^ 

pOO 1 

$D = y (12) 

where the spectral heat flux is given by 

<^0= E [-^[QsG{oj)rf^''{co,K;d) 

This time and denote the transmission co¬ 

efficients of each mode (cj, k) between the source and 
the gate and between the gate and the drain for both 
polarization states j = s, p. In the above relation &ij de¬ 
notes the difference of functions Q{uj^Ti) and Q{uj,Tj). 
According to the N-body near-held heat transfer theory 



FIG. 5 . (a) A classical held-effect transistor is a device with 
three terminals, the source (S), the gate (G), and the drain 
(D). The gate is used to actively control (by applying a voltage 
V bias on it) the apparent electric conductivity of the channel 
between the source and the drain, (b) A radiative thermal 
transistor. A membrane of an IMT material (VO2) acts as 
the gate between two silica (Si02) thermal reservoirs (source 
and drain). The temperature (Tq) of the gate is in steady 
state between the temperatures of the source and drain (Ts 
and Td). and are the radiative heat huxes received by 
the drain and emitted by the source, respectively. 


presented in Ref.^^, the transmission coefficients 

and of the energy carried by each mode written in 

terms of optical reflection coefficients pEj (E = S, D, G) 
and transmission coefficients tej of each basic element 
of the system and in terms of reflection coefficients Pefj 
{E = S,^G and T = S,D, G) of couples of elementary 
elements^ 

T^^^{cO,K,d) 

^ 4 I TG,j 1^ Im(ps,j)Im(pD,j)e-^T'‘^ 

I 1 - PSG,iPD,ie“27<i |2| 1 _ psjpGje-^'f'^ ^ ’ (14) 

t-G/D/, , ^ 4Im(pgG,j)Im(pDj)e-^^‘^ 

/• [LU^rl^a) ^—2'yd 12 

I 1 - PSGjTDje 

introducing the imaginary part of the wavevector normal 
to the surfaces in the multilayer structure 7 = lm{kzo)- 
Similarly the heat flux from the source towards the gate 
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reads 


<^ 3 = E f-$^[QDGit0)T^^^ic0,K;d) 

,=trp}^ ( 27 r) (^ 5 ^ 

+ eGsH7f/®(a;,«;d)] 

where the transmission coefficients are analog to those 
defined in Eq. (iMl and can be obtained making the sub¬ 
stitution S D. At steady state, the net heat flux re¬ 
ceived/emitted by the gate which is just given by the 
heat flux from the source to the gate minus the heat flux 
from the gate to the drain vanishes, i.e. or 


~ = 0 - (lb) 


This relation allows us to identify the gate equilibrium 
temperature (which is not necessary unique because 
of the presence of bistability mechanismsl^ for given 
temperatures Ts and Td- Note that out of steady state, 
the heat flux received/emitted by the gate is = 
~ 7^ 0- If < 0 (^G > 0) an external flux is 

added to (removed from) the gate by heating (cooling). 

To illustrate the different operating modes of the 
NETT we consider now our system of a silica source 
and a silica drain with a VO 2 gate in between. We set 
Ts = 360 K and T]j = 300 K and choose a separation 
distance between the source and the gate and between 
the gate and the drain of d = 100 nm. The thickness of 
the gate layer is set to d = 50 nm. Therefore, the NETT 
operates in the near-held regime (of course, it would also 
work in the far-held regime). The equilibrium tempera¬ 
ture of the gate is obtained by solving the transcendental 


equation (16). Here we hnd = 332 K which is close 


to the critical temperature Tc ~ 340 K of VO 2 which 
means that in the steady-state situation the gate is in its 
insulating phase. In this situation the surface modes of 
the gate and the source and the drain can couple which 
results in large heat huxes 4>s and 4 >d. On the other 
hand, when the temperature of the gate is increased by 
external heating to values larger than Tc then VO 2 un¬ 
dergoes a phase transition towards its metallic phase. In 
this phase we can expect from the discussion of the diode 
that the coupling between the gate and the source and 
the drain will be much less efficient so that 4>s and 4 >d 
can be expected to drop. In the transition regime around 
Tc we model the effective permittivity of VO 2 par¬ 
allel and perpendicular to the optical axis by a simple 
effective medium theory (EMT) ansatz 




_L - f<^T/± + (1 - f)4/± 


(17) 


and by a Bruggemann model (BM) 






+ ( 1 -/) 


,m _ eff 

^ll/-L ^ll/± 


1I/-L 


21 a 


= 0 


(18) 


which is inspired by the measurements and modelling 
in RefP. Here is the permittivity of VO 2 in its 
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FIG. 6. Heat fluxes inside the near-field thermal transistor 
using a VO 2 gate of thickness (5 = 50 nm located at a distance 
d = 100 nm between two massive silica samples maintained at 
temperatures Ts = 360 K and Td = 300 K. By changing the 
flux supplied to the gate the temperature Tq is changed 
which allows to operate the NFTT as a thermal switch, a 
thermal modulator and a thermal amplifier. 


insulating phase and is the permittivity of VO 2 in 
its metallic phase. The filling fraction / is chosen to 
be linear in Tq such that the phase transition starts at 
Tg = 341.2 K and ends at Tg = 345.2 K. Of course, this 
model can be improved by using measurements of the 
permittivities in the transition region. 

Using the above parameters together with the EMT 
model for the permittivity of VO 2 we find the resulting 
heat fluxes 4>s, 4 >d and 4 >g as shown in Eig.[^ These re¬ 
sults illustrate that the NETT can operate as (i) a ther¬ 
mal switch, (ii) a thermal modulator or (hi) a thermal 
amplifier. But let us discuss these functions of the NETT 
in more detail: 

(i) Thermal switching : 

An increase of Tg by about 10 degrees starting from 
Tg = Tq^ leads, as clearly shown in Eig.[^ to a reduction 
of heat fiux received by the drain and lost by the source 
by about one order of magnitude. That means the NETT 
can be used in two operating modes where Tq is slightly 
below or above the critical temperature Tc, where in the 
case Tg < Tc we have a large heat flux to the drain which 
can be assigned as the ”on” mode and for Tg > Tc the 
heat flux towards the drain drops by one order of mag¬ 
nitude which can be assigned as the ’’off’ mode. The 
thermal inertia of the gate as well as its phase transition 
delay defines the timescale at which the switch can oper¬ 
ate. Usually the t hermal inertia limits the speed to some 
microsecond^^^^. 

(ii) Thermal modulation: 

Over a relatively broad temperature range around T^ 
(strip between the vertical lines in Eig.|^ the heat current 
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FIG. 7. Switching mode of radiative transistor: When the 
gate temperature Tq is slightly increased above its critical 
temperature Tc the heat flux received by the drain falls 
from a relatively large value, the so-called ”on” mode, to a 
value which is about one order of magnitude smaller, the so- 
called ”ofF’ mode. 


from or towards the gate remains relatively small 
compared to the and (i-o. ~ Within 

this temperature range the flux received by the drain or 
lost by the source can be modulated. A much larger 
modulation of fluxes can be achieved with the NFTT 
when setting Ts and Td such that « Tc. Then the 
flux can be modulated over one order of magnitude by a 
small temperature change of Tq. 

(hi) Thermal amplification: 

The most important feature for a transistor is its ability 
to amplify the current or electron flux towards the drain. 
In the region of phase transition around Tc we see in 
Fig. that an increase of Tg leads to a drastic reduction 
of flux received by the drain. As described by Kats et al. 
in Ref.SSl this behavior can be associated (in far-field) to 
a reduction of the thermal emission. This corresponds to 
a negative differential thermal conductance (NDTC) as 
described for SiC in Ref.^^and illustrated for our configu¬ 
ration in Fig.[^ Having a NDTC is the key for having an 
amplification which is defined as (see for example Ref.^^ 


d^G 



where 




S/D 


9<Fs/d 
OTg ■ 


(19) 


( 20 ) 


It can be easily seen from Fig. [^that a = 1/2 for Tg 
outside the phase transition region where the material 
properties of VO 2 are more or less independent of Tg, 


FIG. 8. Illustration of the negative differential conductance 
(NDTG) between a VO 2 layer 6 = 50 nm thick and a mas¬ 
sive Si 02 sample separated by a vacuum gap of d = 100 nm. 
The flux $ exchanged between both media is in the phase- 
transition region T 2 G [341.2 K, 345.2 K] a decreasing function 
of the VO 2 temperature T 2 . Outside this transition region 
the flux 4> increases with temperature. Here the temperature 
of Si02 is fixed to 300 K. 
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FIG. 9. Efficiency of the mode coupling in (cj, K:)-space for 
a Si02-V02-Si02 system (d = 50nm and d = 100nm). (a) 
and (b) with VO 2 in its insulating state, (c) 

and (d) with VO 2 in its metallic state. Wien’s 

frequency (where the heat transfer is maximal) at T = 340 K 
is cjwien ~ 1.3 X 10^^ rad/s. 


since in this region. On the other hand, in¬ 

side the transition region of VO 2 that means for temper¬ 
atures around Tc the material properties of VO 2 change 
drastically showing a NDTC which leads to an amplifi¬ 
cation of about 4 regardless of using the EMT or BM for 
modelling the permittivity in the transition regiorP^. 

To explain the strong variations of heat flux observed 
inside our NFTT, let us examine hereafter the transmis¬ 
sion coefficients (only p polarisation) of the energy car- 



































ried by the modes (cj, k,) through such a three-body sys¬ 
tem which are plotted in Fig. When the gate is in 
its insulating state, which represents the exchange 

between the source and the drain mediated by the pres¬ 
ence of the gate [see Fig. M^a)], and which corre¬ 


sponds to the exchange between the couple source-gate 
treated as a unique body and the drain [see Fig. ^b)] 
shows an efficient coupling of modes between the differ¬ 
ent blocks of the system around the resonance frequen¬ 
cies cjsppi ~ 1 X 10 ^^ rad/s and co’spp 2 ~ 2 x 10 ^^ rad/s 
of surface waves (surfaces phonon-polaritons) supported 
by both the source and the drain. Below Tc all parts of 
the system support surface waves in the same frequency 
range close to the thermal peak frequency cjwien ~ 
1.3 X 10^^ rad/s. The anti-crossing curves which appear 
in Fig. I^a) and (b) result from the strong coupling of 
silica surface phonon-polaritons (SPPs) and the surface 
waves (symmetric and antisymmetric ones) supported by 
the thin VO 2 layer. Beyond Tc the gate becomes metal¬ 
lic and it does not sup^rt surface waves anymore. In 
this case 7^^^^ [see Fig. ^c)] vanishes owing to the field 
screening by the gate. Moreover, as is clearly shown in 
Fig. [^d), the coupling of modes between the source-gate 
and the drain at the frequency of surface waves is less 
efficient for the large parallel values of k reducing so the 
transfer of heat towards the drai n, i.e. the number of 
participating modes decrease^^^^^. 


IV. VOLATILE RADIATIVE MEMORY 



FIG. 10. Sketch of a radiative thermal memory. A membrane 
made of an IMT material is placed at a distance d from a 
dielectric layer. The system is surrounded by two thermal 
baths at different temperatures Tl and Tr. The temperature 
T 2 can be increased or reduced either by Joule heating by 
applying a voltage difference through a couple of electrodes 
or by using Peltier elements. 


of simplicity we assume that the separation distance d 
is large enough compared to the thermal wavelengths 
[i.e. d ^ max(ATi = c/i/(/cbT^), i = 1,2,L,R ] so that 
near-field heat exchanges can be neglected (a memory 
using the near-held effect was proposed quite resently in 
Ref.l^. In this case we obtain 


Manipulating radiative heat huxes is only a hrst step 
towards the development of a contactless technology for 
the thermal management of systems at tha nano and 
macroscale. In this last section, we examine the pos¬ 
sibility to store information or energy for arbitrary long 
time using thermal photons by introducing the concept 
of a radiative thermal memory. The concept of thermal 
memory is closely related to the thermal bistability of a 
system that means to the presence of at least two tem¬ 
perature set’s which lead to a steady stat^. 

To illustrate this, let us consider a concrete system as 
depicted in Fig. composed by two parallel homoge¬ 
neous membranes made of VO 2 and Si 02 . These slabs 
have hnite thicknesses and S 2 and are separated by a 
vacuum gap of distance d. The left (right) membrane is 
in contact with a thermal bath having a temperature Tl 
(Tr), where Tl 7 ^ Tr. In a practical point of view, the 
field radiated by these baths can be produced by two ex¬ 
ternal blackbodies. The membranes themselves interact 
on the one hand through the intracavity fields and on the 
other with the thermal baths which can be thought to be 
produced by external media. In that sense, the system 
is driven by many-body interactions in contrast to the 
diode in section II which is only a two-body system. 

The heat flux across any plane z = z parallel to the 
interacting surfaces can again be evalu ated using Ry- 
tov’s fluctuational electro dynamic For the sake 


(p{z) = 2eoC^ y] 

j=s,p 

0 ={ + ,-} 

( 21 ) 


Jo 2 tt J {2TrP 


4>kz0 


UJ 






where the field correlators' 

(w, k) = ^{[Ef (w, K)Ef K)+Ef K)Ef {w, k)]) 


( 22 ) 

of local field amplitudes in polarization j can be ex¬ 
pressed (see the supplemental material iiP^ for details) 
in terms of reflection and transmission operators 91^ and 
of the layer i toward the right (+) and the left (—), 
as 


G:7+=5(Ti)(l-|9i+ P-|T+|2) 

- |2 I rr- |2\ 


c:-’-=5(T2)(i-|9i: 
£+’+ = S{Tl), 

= 5(Tfl), 


■2 I I 


(23) 


with 


5(r)=7r^0(c^,r)Re - 


eoc^ 


(24) 


Here kz denotes the normal component of wave vector 
in the medium of consideration. Using expression (21) 
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we can calculate the net flux = (/:?( 0 ) — [resp. 

$2 = ^{d^ 62 ) — (f{d)] received by the first (resp. second) 
membrane. 



FIG. 11. Trajectories of temperatures (i.e. phase portrait) 
for different initial conditions in the plane (Ti,T 2 ) in a two 
membrane Si02/V02 system with (5i = = 1 /im. The blue 

dashed and red solid lines represent the local equilibrium con¬ 
ditions $1=0 and $2 = 0 of each membrane. The green (red) 
points denote the stable (unstable) global steady-state tern- 
peratures, = (328.03K, 337.77K), = 

(328.06K, 338.51 K), and (Tp\T 2 ^®^) = (324.45K, 341.97K). 
The red arrows represents the vector field $. The tempera¬ 
ture of thermal reservoirs are Tl = 320 K and Tr = 358 K. 

The time evolution of temperatures Ti and T 2 of the 
two membranes are solution of the following nonlinear 
coupled system of differential equations 

= $ + Q (25) 

where we have introduced the vectors T = 
(Ti(t),T 2 (t))\ ^ ^ (^i(Ti,T 2 )//i,$ 2 (Ti,T 2 )// 2 )\ 

and Q = (Qi(^i//i, Q2^2/^2)^- Here Qi {i = 1,2) is 
the power per unit volume which can be added to or 
extracted from both membranes by applying a voltage 
difference through a couple of electrodes as illustrated 
in Fig. or by using Peltier elements. Furthermore, we 
have introduced the thermal inertia of both membranes 
as li = Cipidi^ where Ci and pi are the heat capacity 
and the mass density of each material. By writing down 
this set of equations we have neglected any temperature 
variation inside the membranes which is a very good 
approximation given that the conductivity inside the 
membranes is much larger than between the membranes. 
When assuming that no energy is directly added to or 
removed from the membranes, then Q = 0 . In this case, 
the steady-state solution is given by ^ = 0 . Hence 
$1 and $2 vanish for the same couple of temperatures 
(T^^st“St) ^dst-st)^^ Considering for an instant mem¬ 
brane 2 only, then the existence of two equilibrium 
temperatures where the net flux vanishes ($2 = 0 ) 
implies that $2 niust have a maximum or a minimum 


between these two temperatures. Hence, this requires 
a negative differential conductive behavior for this 
membrane which was shown in the previous section for 
VO 2 . For the whole system it is therefore a precondition 
to have at least one membrane which exihibits negative 
differential conductive behavior in order to have two 
couples Qf steady-state temperatures. 

In Fig.[^ we show the time evolution of the Si 02 -V 02 
system without external excitation (i.e. Q = 0 ) when 
the thermal inertia Ii of both membranes are compara¬ 
ble (i.e. ~ ^ 2 ) or very different (i.e. ^ 62 ). The 

trajectories (the thick pink and turquois lines) are ob¬ 
tained by solving Eq. (26) using a Runge-Kutta method 


with adaptative time steps choosing different initial con¬ 
ditions. In this figure, the dashed blue (solid red) line 
represents the local equilibrium temperatures for the first 
(second) membrane that is the set of temperatures cou¬ 
ples (Ti,T 2 ) which satisfy the condition $i(Ti,T 2 ) = 0 
[$ 2 (^ 1 , T 2 ) = 0]. The intersection of these two lines 
define the global steady-state temperatures of the sys¬ 
tem where $ = 0 . Only two of three equilibrium points 
T = {T^^^T 2 ^Y that appear in Fig. 11 are stable. These 


points can be identified by applying a small perturbation 
6 T = ((5Ti(t), 6 T 2 {t)y on them and look at the relaxation 
process. Then, the time evolution of each perturbation 
is driven by the following equation 


dt{ST)=J^ST, 


(26) 


where 


J = 


1 d^i 1 d^i 
h dTi h dT2 
1 ^*^2 1 ^*^2 
I2 dTi I2 dT2 


(27) 


is the Jacobian matrix associated to the dynamical sys¬ 
tem. Its solution reads 


6 T{t) = JT(0)exp(Jt). (28) 

Hence, the eigenvalues of J allow to determine the 
stability of thermal state. In a two membrane Si 02 -V 02 
system with Si = 62 = 1 /im and with reservoirs temper¬ 
atures Ti^ = 320 K and Tr = 358 K the stability analysis 
leads to 



TP{K) 

Ai 

A 2 

stability 

328.03 

337.77 

-2.58 

-1.37 

stable 

328.06 

338.51 

-2.51 

3.69 

unstable 

324.45 

341.97 

-2.31 

-4.70 

stable 


where we have used the volumetric densities and the heat 
capacities pi = 2180kgm“^, p 2 = 4339kgm“^, Ci = 
750Jkg“^K-i and C 2 = 690Jkg“^K-i. 

From this result, it becomes apparent that it is possi¬ 
ble to use the Si 02 -V 02 system as a thermal memor^ffH. 
Indeed, the two stable equilibrium temperatures can be 
viewed as two thermal states ” 0 ” and ” 1 ” of the system. 
To switch from one thermal state to the other, we need 
to add or extract power from the system. In the follow¬ 
ing we describe this writing-reading procedure. To this 
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Time (s) 


FIG. 12. (a) Hystereis of the VO 2 membrane temperature 
during a transition between the thermal states ”0” and ” 1” in¬ 
side a two membrane Si 02 /V 02 system with (5i = = 1 /xm. 

The volumic powers supplied and extracted from the VO 2 
layer during a time interval Ati = 0.4 s and At 2 = 1.5 s are 
Q 2 = 10“^Wmm“^ and Q 2 = —2.5 x 10“^Wmm“^, respec¬ 
tively. The writing time of state ”1” (”0”) from the state 
”0” (”1”) is At(0 ^ 1) = 4s (At(l ^ 0) = 8s). (b) Time 
evolution of temperature Ti(t) and T 2 (t) of Si02 and VO 2 
membranes. The thermal states ”0” and ”1” can be main¬ 
tained for arbitrary long time provided that the thermostats 
(Tl = 320 K and Tr = 358 K) remain switched on. 


end, we consider the Si 02 -V 02 system made with mem¬ 
branes of equal thicknesses 5i = 82 = 1 /xm which are 
coupled to two reservoirs of temperatures Tp = 320 K 
and Tr = 358 K. Let us define ”0” as the thermal state 
at the temperature T 2 = min(T 2 ^\ T 2 ^^). To make the 
transition towards the thermal state ”1” the VO 2 mem¬ 
brane must be heated. 

Step 1 (transition from the state ”0” to the state ”1”): 
A volumic power Q 2 = 10“^Wmm“^ is added to this 
membrane during a time interval Ati ~ 0.4 s to reach a 
region in the plane (Ti,T 2 ) [see Fig.j^a)] where all tra¬ 
jectories converge naturally (i.e. for Q 2 = 0) after some 
time toward the state ”1”, the overall transition time is 
At(0 1) = 4s Fig. [T^b)]. 

Step 2 (maintaining the stored thermal information): 
Since the state ”1” is a fixed point, the thermal data 
can be maintained for arbitrary long time provided that 
the thermal reservoirs are switched on. This corresponds 
basically to the concept of volatile memory in electronics. 


Step 3 (transition from the state ”1” to the state ”0”): 
Finally, a volumic power Q 2 = —2.5 x 10“^ Wmm“^ is 
extracted from the VO 2 membrane during a time inter¬ 
val At 2 ~ 1.5 s to reach a region [below T 2 = 338 K in 
Fig. p!^a)] of natural convergence to the state ”0” . In 
this case the transition time becomes At(l ^0) = 8 s. 
Compared with its heating, the cooling of VO 2 does not 
follow the same trajectory [see Fig. [T^a)] outlining the 
hysteresis of the system which accompanies its bistable 
behavior. To read out the thermal state of the system 
a classical electronic thermometer based on the thermo 
dependance of the electric resistivity of membranes can 
be used. 


V. CONCLUSION 

In this paper we have summarized the very recent ad¬ 
vances made in optics to manipulate and store the ther¬ 
mal energy with contactless devices. We have demon¬ 
strated the feasability for contactless thermal analogs 
of diodes, transistors and volatile memories. These de¬ 
vices allow for contactless management of heat flows at 
macroscale and at subwavelength scales. These results 
pave the way for a novel technology of thermal manage¬ 
ment. It also suggests the possibility to develop contact¬ 
less thermal analogs of electronic devices such as ther¬ 
mal logic gates for processing information by utilizing 
thermal photons rather than electrons. They also could 
find broad applications in MEMS/NEMS technologies, 
to generate mechanical work by using microresonators 
coupled to a transistor as well as in energy storage tech¬ 
nology to store and release heat upon request. 
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